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ABSTRACT 

Aims. Grain growth has been suggested as one possible explanation for the diminished dust optical depths in the inner regions of 
protoplanetary "transition" disks. In this work, we directly test this hypothesis in the context of current models of grain growth and 
transport. 

Methods. A set of dust evolution models with different disk shapes, masses, turbulence parameters, and drift efficiencies is combined 
with radiative transfer calculations in order to derive theoretical spectral energy distributions (SEDs) and images. 
Results. We find that grain growth and transport effects can indeed produce dips in the infrared SED, as typically found in observations 
of transition disks. Our models achieve the necessary reduction of mass in small dust by producing larger grains, yet not large enough 
to be fragmenting efficiently. However, this population of large grains is still detectable at millimeter wavelengths. Even if perfect 
sticking is assumed and radial drift is neglected, a large population of dust grains is left behind because the time scales on which they 
are swept up by the larger grains are too long. This mechanism thus fails to reproduce the large emission cavities observed in recent 
millimeter-wave interferometric images of accreting transition disks. 
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1. Introduction 

The evolution of circumstellar disks, the birthplaces of planets, 
is still enigmatic, even though the pioneering theoretical work 
on this topic was started almost 40 years ago by Lynden-Bell 
& Pringle (1974). Viscous and/or gravitational stresses are the 
drivers of the disk accretion flows, which can be traced indi- 
rectly by the observed disk lifetimes and accretion rates (e.g., 
Hartmann et al. 1998; Sicilia-Aguilar et al. 2006; Hernandez 
et al. 2007; Fedele et al. 2010). 

While the general trends of declining accretion rates and disk 
masses can be explained by viscous evolution, a sub-set of ob- 
jects called transition disks remain mysterious. These objects 
appear dust-depleted in their inner regions, while the outer re- 
gions resemble normal circumstellar disks (Strom et al. 1989; 
Skrutskie et al. 1990; Calvet et al. 2002; Espaillat et al. 2007, 
2010; Andrews et al. 201 1). The sizes of the dust cavities range 
from a few to more than 70 AU (e.g., Pietu et al. 2006; Hughes 
et al. 2007; Brown et al. 2008; Hughes et al. 2009; Andrews et al. 
2009; Brown et al. 2009; Isella et al. 2010a; Andrews et al. 2010; 
Isella et al. 2010b; Andrews et al. 2011). The gas content of the 
cavities is still largely unknown. It may also be reduced com- 
pared to the outer disk regions as found for example by Najita 
et al. (2010), Dutrey et al. (2008), orLyo et al. (2011); yet other 
works such as Pontoppidan et al. (2008) or Salyk et al. (2011) 
do detect gas inside the dust cavities. 

The time scale of this transition phase of disk evolution is 
estimated to be of the order of a few times 10 5 years (Skrutskie 
et al. 1990; Hartigan et al. 1990). This estimate, however, is 
based on the transition disk frequency, which is still only a lower 
limit due to the lack of spatial resolution. Classifying transition 
disks only based on the SED can be misleading because steep 
decreases in the dust surface density can easily be missed due to 



the presence of small dust particles. Some mechanisms proposed 
to explain observations of transition disks include planet-disk 
interactions (e.g., Rice et al. 2003; Zhu et al. 2011), or photo- 
evaporation (Clarke et al. 2001; Alexander et al. 2006; Ercolano 
et al. 2008). Recent imaging of disks by Andrews et al. (201 1) re- 
vealed a higher fraction (> 20%) of disks with large cavities for 
the mm-brightest sources. Higher fractions of transition disks, 
indicating longer disk clearing time scales, are more difficult to 
explain by photoevaporation. However, the dust emission signa- 
ture of photoevaporating disks has not yet been self-consistently 
modeled, treating grain growth physics and dust-gas feedback. 
It is therefore unclear whether photoevaporating disks leave a 
dust rich or a completely dust and gas free cavity behind (e.g., 
Alexander & Armitage 2007; Garaud 2007). 

The formation of large disk cavities by viscous evolution is 
generally problematic, even with the assumption that photoevap- 
oration or some other mechanism is able to decouple the outer 
and inner disc at larger radii. Indeed the viscous time at 35 AU 
for typical values for the viscosity parameter 1 is 3.4 Myr. On the 
other hand, even at 35 AU, the time scales for grain growth and 
radial drift are only a couple of thousand years. Any mechanism 
that triggers significant changes in the dust evolution (e.g., the 
emergence of a pressure maximum, see Pinilla et al. 2012) could 
therefore quickly induce observational signatures. The fact that 
grains grow and consequently become more mobile due to radial 
drift is likely an important part of the solution to this problem. 

Recent observations by Kraus & Ireland (2012) found a pos- 
sible planetary-mass companion inside the cavity of the LkCa 
15 transition disk. It remains to be shown whether planets cause 
gaps/pressure bumps or vice versa (e.g., Kretke & Lin 2007; 



1 Assuming or, = 3 x 10 3 . Both observational and theoretical works 
indicate values between 10~ 3 and 10~ 2 . 
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Brauer et al. 2008; Pinilla et al., in prep.). Determining the ori- 
gin of transition disks is therefore one of the most fundamental 
issues in our efforts to forge a better understanding of planet for- 
mation. 

It has been suggested, although not demonstrated in any de- 
tail, that the growth and radial transport of dust could potentially 
explain the observed signatures of transition disks with large in- 
ner holes (e.g., Dullemond & Dominik 2005; Tanaka et al. 2005; 
Najita et al. 2008; Pontoppidan et al. 2008). The observations in- 
dicate a decreased optical depth in the inner regions of transition 
disks both in the IR as well as in mm-observations, which does 
not necessarily imply diminished dust densities. Grain growth 
therefore offers two (related) pathways towards this end. First, 
particles decouple from the gas as they grow, which causes them 
to spiral inwards (Weidenschilling 1977; Nakagawa et al. 1986). 
This way, the dust optical depths are reduced by the actual re- 
moval of dust mass. Second, grain growth itself causes a de- 
crease in the dust opacities, since larger particles emit less ef- 
ficiently, which also decreses the otpical depth. Dust grains of 
sizes beyond a few centimeters become basically invisible to 
observations. However, the edges of transition disks seem to 
be relatively abrupt; they have been modeled as step functions 
(Andrews et al. 201 1; Isella et al. 2012) or steep power-law pro- 
files Isella et al. (2010b, 2012). This suggests that the environ- 
ment for grain evolution must abruptly change at this point in 
the disk. Pressure bumps or even gaps opened by planets might 
be possible explanations (Lin & Papaloizou 1986; Zhu et al. 
2011). Another possibility could be the outer edges of dead- 
zones (Gammie 1996). Dead zones are regions where the ioniza- 
tion fraction of the disk drops below the critical value needed to 
drive the magneto-rotational instability (MRI, Balbus & Hawley 
1991), the widely accepted source of turbulent viscosity. 

In this work, we want to investigate whether grain growth 
and transport - either alone or aided by a dead zone - can be 
the cause of the observed transition disk signatures. In Section 2 
we will discuss some of the equations which are crucial for 
the understanding of our modeling results. Section 3 describes 
the model setups and assumptions. The resulting simulation 
outcomes and simulated observations are shown and explained 
in Section 4. Our findings are discussed and summarized in 
Section 5. 



2. Background 

In the following section, we will summarize some of the results 
of Birnstiel et al. (2012, hereafter BKE12), which motivated the 
simulations presented in this paper and will help to understand 
the results. It was found in BKE12 that the upper end of the grain 
size distribution can be limited by two effects, namely fragmen- 
tation and radial drift. In the former case, due to the fact that 
typical impact velocities increase with grain size, the grains can 
only grow until they reach a size at which fragmentation sets in. 
The size limit for turbulent relative velocities was found to be 

E m 2 

a fmg * 0.08 (1) 

Here, Uf denotes the fragmentation threshold velocity and p s the 
specific density of the dust grains. E g is the gas surface den- 
sity, c s the sound speed and a t the turbulent viscosity parameter 
(Shakura & Sunyaev 1973). The other size limit is due to ra- 
dial drift: grains can only exist at a given radius if growth from 
smaller sizes resupplies them as fast as radial drift removes them. 



In this case, the upper end of the size distribution can be approx- 
imated by (see BKE12) 

a dnft *0.35^(^V, (2) 

where Ed is the dust surface density, the Keplerian frequency 
and y the absolute value of the power law index of the gas pres- 
sure P, 



and we have introduced a parameter for the efficiency of radial 
drift, Ed where E& = 1 corresponds to the fiducial literature 
value. In the limit flf lag < Adrift, fragmentation of the dust grains 
is the relevant size limit. In this case, small dust is constantly 
resupplied in the form of fragments. In the opposite case, when 
the radial drift barrier is the growth limiting factor, grains are 
removed by radial drift even before they reach sizes at which 
fragmentation starts to become important. 

From Eqns. 1 and 2, it can be seen that the quantities which 
most clearly determine the time evolution of the size limits are 
Eg and Ed, respectively, because the other quantities are not ex- 
pected to change by orders of magnitude with time at a given 
radius. However, it should be noted that in this study, we do not 
track the viscous evolution of the gas surface density. This would 
lead to a further reduction of the maximum grain size. A fixed 
gas surface density is therefore in tune with having a best-case 
scenario to test the potential of the proposed mechanism. The 
fact that radial drift causes the dust surface density to decrease 
faster than the gas surface density, even if viscous evolution were 
to be included, means that the drift limit can become more im- 
portant as the disk evolves. 

The dependence of Eqns. 1 and 2 on E e and Ed also causes 
both size limits to decrease with radius, as observed by Banzatti 
et al. (2011) and Guilloteau et al. (2011). If the largest grain 
size at a given radius is limited by radial drift, dust collision 
velocities are not high enough to cause fragmentation (BKE12). 
The drifting grains therefore sweep up the small dust without 
resupplying it by fragmentation. This means that drift-limited 
grain growth could naturally lead to conditions where the inner 
regions are devoid of small dust (due to lack of fragmentation), 
while small grains are still present in the outer disk (due to lower 
grain size limits in the outer regions). 

The most important and uncertain parameters are Uf, the col- 
lision velocity at which fragmentation sets in, and the turbu- 
lent state of the disk, which is described by the parameter a t . 
Estimates for Uf range from a few up to 35 m s _1 (see, Wada 
et al. 2008; Paszun & Dominik 2009; Blum & Wurm 2008, and 
references therein). We use values of 3 and 10 m s _I , because 
icy aggregates are expected to fragment only at larger velocities 
compared to the ~lms _l , found for silicate dust grains (Blum 
& Muench 1993; Blum & Wurm 2008). 

3. Modeling approach 

3.1. Dust evolution 

In this work, we use a vertically averaged dust evolution code 
which tracks the radial and size evolution of the dust surface den- 
sity of the disk. Effects of coagulation, fragmentation, and cra- 
tering as well as radial drift, gas drag, and turbulent mixing are 
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Fig. 1. Initial conditions of the simulations without a dead zone. 
Solid and dashed lines denote the gas and dust surface density in 
g cm' 2 , respectively. The dash-dotted line denotes the mid-plane 
temperature. 

taken into account. For any details of the dust evolution model, 
we refer to Birnstiel et al. (2010). 

We model both active disks, where the complete disk is MRI 
active as well as disks with MRI-inactive regions, the so-called 
dead zones. The properties and even the mere existence of dead 
zones is still an active matter of debate. In order to investigate 
whether a dead-zone can produce the observational appearance 
of a protoplanetary disk, we have to make some simplifications 
and assumptions. Therefore, we include the two most prominent 
effects of a dead zone which can alter the evolution of dust, 
namely the decreased amount of turbulence and increased gas 
surface density in the MRI inactive regions. 

For the gas disk, we use a fixed surface density profile with 
a constant mass accretion rate up to a characteristic radius r c , 



2 g (r)oc 



Q k 




• exp 



(4) 



a t c; \ v r 

where r; n is the innermost radius of the disk. In order to mimic 
the effects of a dead zone, we use an alpha profile of 



a t {r) = or D - (a D - a A ) ■ 



: exp 



1 



v Ar 

■ exp ( 



Ar 



if r < r dz 
if r > r dz 



(5) 



which represents a smooth transition at r& z from the turbulence 
parameter in the dead zone au to the active one a a over a transi- 
tion width Ar, which we arbitrarily set to 1 AU. This jump in a t 
causes also a jump in the gas surface density profile, as can be 
seen from Eq. 4. 

The initial condition for the dust surface density is given by 
a constant dust-to-gas ratio in jjm sized grains. The temperature 
profile is derived from two-dimensional radiative transfer calcu- 
lations (see Sect. 3.3) and is kept fixed throughout the simula- 
tion. As an example, the initial gas and dust surface densities 
and the initial temperature profile of simulation A2 is shown in 
Fig. 1. 




Fig. 2. Vertically integrated dust surface density distribution af- 
ter 5 Myrs of evolution for the simulations without a dead zone, 
A2 (top) to A5 (bottom). The parameters for the simulations are 
shown in Table 1 . The solid red line denotes the growth barrier 
set by grain fragmentation, the dashed red line the size limit due 
to radial drift. 



3.2. Vertical structure 

Following the results from Dubrulle et al. (1995), the grain evo- 
lution code used in this work assumes a Gaussian vertical dis- 
tribution for the dust, where the scale heights depend on the 



grain size (see Birnstiel et al. 2010). This assumption is very 
accurate at the mid-plane, but it deviates in the surface layers 
of the disk. It can therefore be used for the dust size evolution, 
which mostly depends on the mid-plane values. However, to de- 
rive proper observables from the simulation outcome, the two- 
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between vertical mixing and dust settling 

<5 I Pa 



dpd 
dt 



d_ 
'dz 



Pd«sett - D d p g 



Fig. 3. Vertically integrated dust surface density distribution af- 
ter 5 Myrs of evolution for simulations D23J105 (top), D25_M85 
(second panel), LDE (third panel) and COAG (bottom panel). The 
corresponding parameters are shown in Table 1. The solid red 
line denotes the growth barrier set by grain fragmentation, the 
dashed red line the size limit due to radial drift. 



dimensional dust density distribution needs to be reconstructed 
from the dust surface densities. 

To this end, the vertical structure is assumed to be station- 
ary and isothermal and we solve numerically for the equilibrium 



VPg 

where pa and p g are the dust and gas densities. 
D A 



(6) 



1 + St 2 



(7) 



is the dust diffusivity (Youdin & Lithwick 2007); we assume that 
the gas diffusivity equals the gas viscosity v g . St is the particle 
Stokes number and 



3 ■sjnmzti^ 
8 V2pg n a 2 c s 



(8) 



is the settling velocity for spherical grains in the Epstein regime 
(Nakagawa et al. 1986). We denote the grain masses and radii 
as m and a, respectively. The resulting two-dimensional density 
distribution of each grain size is then used in the radiative trans- 
fer calculations. 



3.3. Radiative transfer calculations 

The dust evolution models described above need to be con- 
verted into synthetic data products to assess how well they can 
reproduce the key observational properties of transition disks. 
For each model, we compute a broadband SED and high res- 
olution millimeter-wave continuum image following the basic 
procedures outlined by (Andrews et al. 2009, 2011). The dust 
density structure pj for each grain size a was re-sampled onto 
a fixed grid in spherical coordinates, with high resolution re- 
finements near the disk midplane and around regions of strong 
optical depth gradients (i.e., at the inner edge near 0.1 AU and 
the outer boundary of the dead zone). Absorption and scattering 
opacities for each grain size were calculated with a Mie code, 
assuming a population of segregated spheres with the optical 
constants and material compositions advocated by Pollack et al. 
(1994) (see also Andrews et al. 2012). We assume that stellar 
irradiation is the only relevant heating source, and set the stellar 
photosphere to have properties that are representative of typical 
transition disk hosts: T e s = 4300 K, R* = 2.5 R© (implying 
L* * 2 L Q ), and M* = 1 M Q . 

The two-dimensional Monte Carlo radiative transfer code 
RADMC (Dullemond & Dominik 2004) was used to simulate the 
propagation of radiation through each model structure and com- 
pute an internally-consistent temperature structure. That process 
was iterated with the dust evolution code for each model so ad- 
justments could be made to the vertical distribution of particles 
of a given size. We adopted a simple convergence criterion such 
that the midplane temperatures did not change more than ~10% 
between iterations (in practice, this amounted to only 1 or 2 iter- 
ations). A ray tracing algorithm was then used on the final model 
structure to construct a theoretical SED and 880 pm continuum 
image for a given disk viewing geometry and distance. For sim- 
plicity, we fixed representative values for the latter; d = 140 pc, 
a disk inclination i = 35°, and major axis position angle of 155° 
(measured east of north). The 880 pm model image is used to 
synthesize a simulated interferometric image like those observed 
by the Submillimeter Array (for examples, see Andrews et al. 
2011). 
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3.4. Shortcomings of the modeling approach 

A vertically integrated coagulation code combined with a ver- 
tical settling-mixing equilibrium calculation is a good approxi- 
mation of the vertical dust structure, as long as the disk is not 
layered, as it is in the presence of a dead zone. For example, a 
small amount of the mid-plane dust could reach the surface lay- 
ers and be fragmented due to the high degree of turbulence in 
the active regions. Furthermore, turbulent transport in these re- 
gions can be much stronger than in the mid-plane (Turner et al. 
2010). Due to the fact that our model only uses the (possibly 
very low) mid-plane turbulence parameter, this model represents 
a 'best case' scenario in favor of the removal of small grains. For 
simplicity, the evolution of the gas surface density and also the 
dependence of the MRI turbulence on the dust distribution (e.g., 
Sano et al. 2000) are not treated in this work. 

It is important to note that other obstacles for particle growth 
at smaller sizes have been suggested, e.g., bouncing (Guttler 
et al. 2010; Zsom et al. 2010) or charging effects (Okuzumi 
2009). Observations, however, do find ample evidence of the 
presence of both large (Testi et al. 2001; Natta et al. 2004; 
Rodmann et al. 2006; Ricci et al. 2010) and small grains (e.g., 
Watson et al. 2007, and references therein) in protoplanetary 
disks, which suggest that these barriers have been overcome. The 
question remains whether the radial drift and the fragmentation 
barrier are indeed the upper limits of the size distribution. At 
least fragmentation seems necessary in order to explain the ob- 
served presence of small dust (see Dullemond & Dominik 2005), 
however the expected rate of radial drift is too high to account for 
the observed disk lifetimes (Weidenschilling 1977; Brauer et al. 
2007). On the other hand, different sizes as well as different pro- 
files of the dust and the gas disk (Panic et al. 2009; Andrews et al. 
2012) are indicating that to some extent radial drift is at work in 
circumstellar disks. Possibly reduced drift rates (e.g., Johansen 
et al. 2006) or pressure traps (Pinilla et al. 2012) are needed to 
reconcile theory and observations. To account for these open is- 
sues, we will also present models with reduced or without radial 
drift and without fragmentation (see Sect. 4.1). 

4. Results 

In this section we will discuss the results of two parameter stud- 
ies: one of entirely MRI-active disks with varied levels of tur- 
bulence, and another which mimics MRI-inactive regions (dead 
zones) in which we vary the turbulence in both the MRI-active 
and dead regions as well as the efficiency of radial drift. A sum- 
mary of the models and parameters is given in Table 1 . 

4. 1 . Growth and transport in smooth disk profiles 
4.1 .1 . Evolution of the dust distribution 

In this first parameter study, we investigate how the possible 
transition from a fragmentation dominated («f < aj) to a drift 
dominated (aj < Of) size distribution is reflected in the spectral 
energy distribution (SED) of a typical disk profile. To this end we 
set up different simulations with the same initial conditions and 
vary the constant turbulence parameter a t from 10~ 2 to 10~ 5 in 
simulations A2 to A5, respectively. According to Eq. 1, this vari- 
ation in a t shifts the fragmentation barrier with respect to the 
drift barrier (Eq. 2) by three orders of magnitude in size. Thus 
the maximum size of the dust grains in the simulation with the 
highest a t value of 10 -2 (i.e., simulation A2) is set by fragmen- 
tation, while the drift limit becomes more important as we go to 



lower turbulence values. In this parameter study, we vary only 
the turbulence parameter, however it is important to note that 
the same or similar effects can be achieved by varying the initial 
dust-to-gas mass ratio Ed/Eg (for Ed « E g ), or the fragmentation 
velocity Uf, since these parameters determine the initial ratio of 
the size limits (Eqns. 1 and 2). Changing these parameters, a 
drift-limited size distribution can be achieved, even for a disk 
with relatively strong turbulence. 

Figure 2 shows <x(r, a), the dust surface density distributions 
as function of grain size and radius for the active disk simula- 
tions, A2 through A5, after 5 Myrs of evolution. The quantity 
cr(r, a) is defined as 

rOO 

cr(r, a) = n(r,z,a) madz, (9) 

J — 00 

where m = 4n/3p s a 3 is the particle mass, and n(r,z,a) is the 
dust number density distribution as a function of radius r, height 
above the mid-plane z, and grain size a. Integration of <r(r,a) 
over ln(fl) gives the dust surface density Sd(r). The solid red line 
in Fig. 2 denotes the fragmentation limit (cf. Eq. 1), the dashed 
red line the drift size limit (cf. Eq. 2). 

Going from simulation A2 to A5 (see Fig. 2), it is evi- 
dent that the small grain population becomes increasingly de- 
pleted. The reason for this behaviour is the following: radial 
drift causes the dust surface density to be decreased on shorter 
time scales than the gas surface density. The dust-to-gas ratio 
is therefore decreasing with time. The ratio of the two growth 
barriers, which in this study is initially set by the choice of a t , 
therefore shifts with time from a fragmentation limited distribu- 
tion towards a drift limited one. As described in BKE12, for a 
drift limited size distribution, dust collision velocities induced 
by radial drift decrease with the dust-to-gas ratio. Therefore, at 
later times, the low collision velocities prevent grain fragmen- 
tation. Consequently, small dust particles are not reproduced by 
fragmentation anymore, but are rather swept up by the inward- 
drifting population of larger grains. 

The time scale of this sweep-up can be estimated in the con- 
text of the drift-dominated two-population model of BKE12 (see 
Appendix A) as 

where e is the dust-to-gas mass ratio, td and tj are the sweep-up 
time scales for drift and turbulence induced collision velocities. 
The dependence on temperature T and y is not relevant, con- 
sidering that these values do not change drastically with time 
or radius and the dependence on T is weak. The sweep-up time 
scale is therefore set by the quantities r, e, and a t . The late stages 
of disk evolution after 5 Myrs in Fig. 2 are typically paired with 
dust-to-gas mass ratios of 10~ 4 to 10~ 5 , quite irrespective of 
at- 2 The radial dependence and time evolution of the dust to gas 
ratio is very similar to the results shown in Fig. 8 of BKE12. 

2 As long as the dust distribution becomes drift dominated, this leads 
to a self-regulating process: radial drift removes dust mass, thus lower- 
ing the dust-to-gas ratio. As the dust-to-gas ratio drops, the growth time 
scale and therefore also the drift time scale increases. Hence, the evo- 
lutionary time scale of the dust-to-gas ratio depends on the dust-to-gas 
ratio itself, and to reach a dust evolution time scale of a few Myrs, the 
dust-to-gas ratio needs to be in the range of 10~ 5 to some 10~ 4 . 
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Table 1. Parameters and initial conditions of the simulations 
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Notes. Models A2 through A5 are the completely active disks, models starting with D are the disks with dead zones, LDE stands for low drift 
efficiency, and COAG is a model with neither radial drift nor fragmentation. Other relevant initial conditions are the stellar mass M t = 1 Mq, stellar 
temperature T„ = 4300 K, and/?* = 2.5 Rq. 
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Fig. 4. Spectral energy distributions (top row), normalized radial surface brightness profiles at 880 fim (middle row), and synthesized 
880 fim images (bottom row) corresponding to the Submillimeter Array setup of Andrews et al. (201 1) for the active disk models. 
The line colors correspond to models A2 to A5 and model COAG, as labeled in the left panels. The dashed grey line in the SEDs 
represents the stellar photosphere. The columns from left to right correspond to 0, 1, 2, and 5 Myrs of evolution. The shaded areas 
in the central row depict the resolution limits of current (pre- ALMA) interferometers. Only the results for model A4 are shown in 
the synthesized images, where the contours correspond to 10% levels of the peak brightness and the bar and dashed oval represent 
the size scale and the beam size, respectively. 
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Fig. 5. Same as Fig. 4 but for the models which include dead-zones, corresponding to the simulations labeled in the left panels. 
The dashed vertical line in the central row marks the dead zone radius r^. The synthesized images in the bottom row correspond to 
modelD25_M<95. 



Therefore, the sweep-up time scales approach ~Myrs values 
beyond about 10 AU. This explains why dust at these radii is 
not swept up quickly. The linear dependence with r means that 
the clearing of dust in this picture proceeds from inside-out, as 
typically associated with other disc-clearing mechanisms (e.g., 
Ercolano et al. 2011). In the innermost regions, the drift size 
limit approaches the fragmentation limit. This causes some frag- 
mentation of the largest particles, which is the source of the 
small dust population inside of ~0.5 AU for simulation A5 and 
inside of ~ 1 AU for simulation A4. 



4.1.2. Evolution of the SED 

The resulting SEDs and 880 fim intensity profiles for all the 
MRI-active disk simulations are shown in Fig. 4, where the pan- 
els from left to right show the time evolution of the different 
models. It can be seen that the mid-IR emission of A2 and A3 are 
comparable up to about 20 //m, because even at 5 Myrs the frag- 
mentation barrier in both cases is still at similar or even smaller 
grains sizes than the drift limit. Going to the lower a t values of 



A4 and A5 causes an ever stronger reduction of the mid-IR emis- 
sion, due to the fact that the drift limit becomes significantly 
smaller than the fragmentation limit and the small dust grains in 
the hot inner regions are swept up by the drifting grains. 



The outer regions show a similar effect: the emission beyond 
about 30 /mi stays high in the case of A2 due to the fact that 
the strong level of turbulence prevents the formation of larger 
grains. The dust grains in this case are so small that they are 
well coupled to the gas. Hence, the dust is not strongly drifting, 
which means that more mass is retained in the outer regions and 
the emission therefore stays at a high level. Reducing the turbu- 
lence strength in the fragmentation limited case leads to larger 
particles which start to attain significant inward drift velocities. 
The far-IR flux is therefore decreased, because less dust mass is 
retained in the outer regions of the disk. This explains the strong 
difference between the SEDs of model A2 compared to the other 
models. 



7 



T. Birnstiel et al.: Can grain growth explain transition disks? 



4.1.3. Evolution of the (sub-)mm brightness 

The shape of the final radial intensity profiles in the 10 - 100 AU 
range can be understood by approximating I v oc T 2^ k v . For sim- 
ulations A3, A4, and A5, the largest grain size between about 
7 and 50 AU is given by the drift-limit Eq. 2 and at the same 
time is larger than 880 [im./(2n). This means that the opac- 
ity #880 is roughly k%%q oc a~ oc T y r/Ej. The intensity then be- 
comes I v a:T 2 yr, which for a temperature profile T ~ r~ l l 2 
leads to a flat intensity profile. The entire disk of model A2 
is fragmentation-dominated. This gives rise to a steeper inten- 
sity profile because the presence of large amounts of small dust 
means the opacity does not depend as strongly on the largest 
grain size. For a typical size distribution in the fragmentation 
case, we get 7 v oc ^/E g /r which explains the slope of model A2. 

4.1 .4. Conclusions for the smooth disk models 

Even though the SED signatures of the low turbulence simula- 
tions do resemble the ones of typical transition disks, there is an 
important drawback of this mechanism when it comes to disks 
with resolved holes in the (sub-)millimeter wavelength range: as 
pioneered by Brown et al. (2009) and confirmed by other stud- 
ies such as Andrews et al. (2011), many of the transition disks 
(but also some of the disks without typical transition disk-like 
SEDs) were found to have inner dust cavities with sizes up to 
70 AU. While it is easy to "hide" such cavities in the SED by rel- 
atively small amounts of small dust, the resolved optically thin 
(sub-)millimeter observations unambiguously reveal these fea- 
tures, as shown e.g., by Isella et al. (2010b) and Andrews et al. 
(2011). 

The mechanism presented here is able to remove the small 
dust population; however, it does so by locking it up in a non- 
fragmenting population of larger dust particles. This population 
of grains contains still enough particles of ^centimeter sizes 
such that they clearly show up in the (sub-)millimeter images. 
So even in the best case (model A4), no cavity is seen in the 
simulated images (Fig. 4, bottom row). The size of the largest 
particles scale with the gas surface density; however, we chose 
a very high disk mass (20% of the mass of the central star) in 
order to reach the largest possible sizes for the given disk model. 
Even in this case, the images shown in Fig. 4 do not resemble 
the observational data of e.g., Brown et al. (2009). In the follow- 
ing section, we will investigate how far this mechanism can be 
pushed by including effects which might be caused by a layered 
structure of the disk. 

4.2. Beyond standard dust modeling 

4.2.1 . Including effects of a layered disk structure 

So far, we have shown that the mechanism presented in this pa- 
per is able to strongly deplete the inner regions of disks in small 
dust grains and that a mid-IR dip in the SED, as found by obser- 
vations (e.g. Skrutskie et al. 1990; Calvet et al. 2002; Espaillat 
et al. 2007, 2010) can be produced. In order to resemble the 
resolved (sub-)mm observations such as Brown et al. (2009), 
Andrews et al. (201 1), or Lyo et al. (201 1), two further features 
are required. First, the dust cavities also need to be observed 
in the (sub-)mm, i.e. for larger grains. The absence of a cavity 
in the simulated images of Fig. 4 indicates that the grain sizes 
found in our simulations are not large enough. Larger surface 
densities or larger fragmentation thresholds can be considered. 
The former option, however, is restricted by the total mass of the 



disk. Second, in our simulations, the dust surface density and the 
largest grain sizes are smooth functions of the radius, giving rise 
to a flat intensity profile (cf. the estimate in Sect. 4.1). This is 
in contrast with the observed dust cavities. These facts motivate 
a model which allows for higher gas surface densities and steep 
changes in the maximum grain size at the cavity edge. 

One candidate for this behavior could be disks with dead 
zones, regions of low ionization which are not MRI-active. 
The inactive regions transport angular momentum less effi- 
ciently (Armitage 2011), which is why they can accumulate 
large amounts of mass, while the outer, active regions still re- 
semble normal disks. This way, the inner, MRI-dead regions can 
have higher surface densities, while the transition from MRI- 
active to MRI-dead gives rise to strong changes in the turbulence 
parameter a t (e.g. Dzyurkevich et al. 2010). 

In this section, we therefore take these two effects into ac- 
count without any detailed modeling of the dead zone. This 
means that neither the vertical structure nor the viscous evolu- 
tion of the dead zone are considered; only the reduced a t and 
the increased gas surface density S g are treated. This is a valid 
approach because we are considering this as a best-case model. 
Other effects such as fragmentation and enhanced radial mix- 
ing in the active surface layers are working against the presented 
mechanism by replenishing small dust grains. 

The situation where a t is higher in the outer parts of the disk 
and lower inside the dead zone can cause the dust distribution to 
be fragmentation limited outside and drift limited inside the dead 
zone. Models D23J105 and D25J105 simulate such a situation 
for rfc = 35 AU, but for different contrasts between the active 
and the dead region. The size of the dead zone is not physically 
motivated, but rather chosen from the typical size of the observed 
transition disk cavities. However, the theoretical models of Bai 
(201 1) can account for dead zones of sizes up to 20 AU. 

A comparison of the SEDs for model D2 3 _M®5 (Fig. 5, purple 
SED) and A2 (Fig. 4, purple SED) highlights how the contrast 
of one order of magnitude in a t changes the long wavelength 
emission in the SED. For our choice of parameters, the outer re- 
gions start to be effected by drift, but fragmentation stays active 
throughout the disk. Only at 5 Myrs, between around 10 AU and 
the outer edge of the dead zone, the drift barrier drops slightly 
below the fragmentation barrier, causing a minor reduction of 
small grains. 

We also compared models with different disk masses (not 
shown), but no significant differences were found for changes 
within a factor of a few. Going to very low disk masses of 
0.005 Mq further reduces the maximum grain size, which is pro- 
ducing relatively more mm- and IR emission in the inner regions. 
Severe changes are found if the contrast in a t between the active 
and the dead regions is increased, because this changes the rel- 
ative importance of the fragmentation and the drift barrier, as 
already discussed in the previous section. The simulated images 
in Fig. 5 (bottom row), however, do not show any signs of a cav- 
ity. The reason for this is that the peak of the size distribution 
reaches sizes of around a centimeter and larger only within a 
few AU (see third plot from top of Fig. 3). Therefore, the surface 
density of dust smaller than ~ 1 cm only decreases smoothly 
inwards of 3-8 AU, depending on the age of the disk. 

In order to confirm this trend, we carried out simulations in 
which the size limits are pushed to even larger values. One rea- 
son for larger grains can be the efficiency of radial drift: simula- 
tions by Johansen et al. (2006) and Bai & Stone (2010) showed 
that the radial drift velocity in the case of a turbulent environ- 
ment can be reduced compared to the laminar, single-size results 
of Weidenschilling (1977) and Nakagawa et al. (1986). In order 
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to artificially change the efficiency of radial drift, we have intro- 
duced the efficiency factor in Eq. 3. From Eq. 2, it follows 
that Adrift then becomes inversely proportional to the drift effi- 
ciency. This means that if radial drift is only Ea = 0.1 times as 
fast as the classical value, flf rag became l/E^ = 10 times larger 3 . 

Our simulation results with reduced drift efficiency (cf. 
model LDE in Table 1) show that particles do reach larger sizes 
in the dead zone (see third panel from top in Fig. 3). Still, the 
grain size is a smooth function of radius inside the dead zone. 
We found that going from 20 to 1 AU, the total surface density 
of grains smaller than one centimeter decreases only smoothly 
by two orders of magnitude. This means that the cavities in the 
mm-images cannot be explained with this mechanism. In addi- 
tion to that, the reduced drift speed also increases the time scale 
on which small dust is swept up (cf. Eq. 1 1). This results in much 
larger amounts of small dust being left over, as can be seen by 
comparing the second and third panel in Fig. 3. It also explains 
the less pronounced infrared dip in the SEDs (cf. Figs 4 and 5). 

The intensity shapes inside the dead zones can be understood 
as in the previous subsection: for all models but D2 3 J105, we re- 
cover a rather flat intensity profile due to the presence of large 
grains and a drift limited distribution. In model D23_M05, the 
effective fragmentation causes a steep increase in the inner re- 
gions, as in model A3. This increase is mostly caused by changes 
in the opacity due to the presence of large amounts of small dust. 

4.2.2. Growth through all barriers 

It is a valid question to ask what disks would look like if our 
current understanding of the growth barriers and radial drift is 
entirely wrong. For this reason we simulated a case where grain 
growth is not inhibited by either fragmentation or radial drift. 
That is, we consider perfect sticking (without any fragmentation) 
and in addition, dust is only allowed to be diffused or dragged 
along by the gas, but not to drift (i.e., E^ = 0). Apart from those 
changes, all other parameters of this model are identical to A55 
(see model COAG in Table 1). The resulting grain size distribution 
is plotted in the bottom panel of Fig. 3 and the resulting SEDs 
and brightness profiles are shown in Fig. 4. 

In Fig. 4, it can be seen that for perfect sticking, the IR dip 
in the SED forms very quickly and is very pronounced, as also 
found by Dullemond & Dominik (2005). While the inner parts 
are cleared of small dust, the growing dust particles are not able 
to sweep up all smaller dust further out in the disk, a population 
of small dust remains, as seen in the bottom panel of Fig. 3. We 
calculated the time scale for the small dust sweep-up in such a 
scenario (see Appendix A), which is proportional to the planetes- 
imal size and to the quantity r/(SdC s ) and for our setup exceeds 
1 Myr beyond 10 AU. This approximation does not include a 
reproduction of small dust, as would be expected for a distri- 
bution of planetesimals, which would further increase the time 
scale. Therefore, even a scenario with perfect sticking and no ra- 
dial drift is not able to reproduce the observed large cavities in 
transition disks. 



3 It should be noted that in the fragmentation dominated limit (af iae < 
Adrift), the drift velocity of the largest particles is decreased if Ea is 
decreased while the particle size does not change. In the drift limit 
(tfdrift < afag), the situation is reversed: the largest particles will grow 
to larger sizes until their drift time scale again equals the growth time 
scale. Thus the particle size is increased while the drift velocity of the 
largest particles is unchanged as long as their Stokes number is below 
unity. 



5. Discussion and conclusions 

In this paper, we have investigated the ability of models of dust 
evolution to explain transition disks. We calculated simulated 
SEDs and 880 fim images from the output produced by a dust 
evolution code for a variety of different initial conditions. For 
clarity, we only considered effects directly induced by growth, 
fragmentation and transport of dust, without pressure traps or 
planet induced gaps. The models were chosen to represent best- 
case scenarios for this mechanism to work, i.e. large gas disk 
masses and no viscous evolution were assumed. The only exten- 
sions to this model we have considered here are the ones which 
influence the grain size limits, i.e., the strength of turbulence and 
the gas surface density. We found that effects of grain growth can 
indeed produce dips in the SED such as found for many transi- 
tion disks (see also Dullemond & Dominik 2005), however they 
fail to reproduce the cavities in the (sub-)millimeter images such 
as found by Brown et al. (2009); Andrews et al. (201 1); Lyo et al. 
(2011). Even in the extreme case where growth proceeds with- 
out any barriers and radial drift is switched off, no large inner 
cavities can be formed within 5 Myrs of evolution. 

This leads to the conclusion that disks with large inner holes 
cannot be caused by grain growth alone. Due to the fact that 
current models of disk photoevaporation (e.g. Owen et al. 2010) 
also fail to explain the inner cavities in the accreting objects, 
we propose that a combination of dust evolution with other ef- 
fects such as pressure bumps, spiral arms or planet induced gaps 
could be the solution of the problem (see Pinilla et al., submit- 
ted to A&A). Important questions to be answered from the ob- 
servational side are: apart from the size of the cavity, is there 
is a distinct difference between disks with small and large cavi- 
ties? What is the gas content inside the cavities? And do different 
sizes of dust show different cavity shapes in the same object (see 
Dong et al. 2012)? Theoretical models will need to investigate 
the trapping mechanisms which effectively shepherds the dust 
outside the cavity, irrespective of the particle sizes. 

Our findings can be summarized as follows: 

- A grain size distribution for which the drift-induced grain 
size limit is smaller than the fragmentation induced size limit 
becomes inefficient in replenishing small dust due to a lack 
of fragmenting collisions. 

- The time scale on which the remaining small dust is swept 
up by the largest, inward drifting grains is proportional to the 
distance to the central star and depends also on the dust-to- 
gas mass ratio and the turbulence strength. The dispersal of 
small dust is therefore from the inside out. 

- Small amounts of small dust can be retained for several mil- 
lion years in the outer parts of the disk (> 10 AU) if the 
dust-to-gas mass ratio is lower than a few times 10~ 4 . 

- Even if the inner regions are cleared of small dust via grain 
growth, most of the dust mass is still detectable by millimeter 
observations. Furthermore, the grain size is a smooth func- 
tion of the stellocentric radius which means that the "observ- 
able" dust surface density decreases only slightly and the 
brightness profile stays approximately flat. This mechanism 
thus fails to reproduce the strong drop in millimeter emission 
which is typically observed in disks with large cavities. 

- The model does not resemble the observations of large cav- 
ities, even if perfect sticking and no radial drift is assumed. 
The observed, sharp division between the inner and the 
outer regions call for a severe change in the disk properties. 
Pressure bumps that are strong enough to decouple both large 
and small dust from the accretion flow are most likely nec- 
essary to explain the observations. 
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Appendix A: Estimation of the sweep-up time 
scales 

To estimate the time scale on which small dust is swept up by 
larger, drifting particles, we consider the two-population model 
of BKE12 according to which the grain size of the drifting par- 
ticles is given by 

fli = — f, (A.l) 

np s yct 

where /d = 0.55 and Vk is the Keplerian velocity. The drift 
velocity of the largest grains is then given by 



(A.2) 



where e is the dust-to-gas mass ratio. Assuming that the large 
particles contain most of the mass of the system (Si ~ Xd), and 
that the small particles are well mixed with the gas (Sti < a t ), 
we can write the change of the surface number density of small 
particles A^o as (see, Birnstiel et al. 2010, Appendix A.2.) 



No 



1 



MoicroiA^oA^i. 



(A.3) 



Approximating the relative velocity m 01 between the two species 
by the drift velocity of the large grains and the collision cross 
section by cr 01 ~ ncP[, we can derive the time scale for the 
sweep-up of small dust by drift-induced collisions 



No 
N 



n ec s y 



(A.4) 



The same calculation but using the turbulent relative velocities 
between the largest and smallest grains 



M 1 



a t a\ p s n 



(A.5) 



according to Ormel & Cuzzi (2007) yields the sweep-up time of 
small particles where collisions are driven by turbulence, 



(A.6) 




10 



T. Birnstiel et al.: Can grain growth explain transition disks? 



Comparing Eqns. A.4 and A. 6 shows that the sweep-up of the 
small dust is driven by turbulent velocities, as long as 

a t >^f d y-£, (A.7) 

which is in our models approximately gives a t > I . 

A similar calculation can be done for the small dust in a dis- 
tribution which is dominated by planetesimals of size ci\. Using 
Eq. A. 3 with the cross section <x 01 ~ jia 2 ^ and the relative veloc- 
ities from the head wind Amoi = yj (2Vk), we derive 

8 V^rps a\r 

Tp = = • (A.8) 

3 At c s y 

For the values of our simulations (model COAG) and a\ = 
10 4 cm, we find that Tp ~ 1 Myrs at 10 AU and increases steeply 
with radius. This agrees well with the fact that the planetesimals 
cannot efficiently sweep up the small grains beyond a few AU, 
as seen in the bottom panel of Fig. 3. 
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